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Abstract 

The quark antiquark potential, and its associated zero point energy, can be extracted from lattice 
measurements of the Wilson loop. We discuss a unique prescription to renormalize the Wilson loop, 
for which the perturbative contribution to the zero point energy vanishes identically. A zero point 
energy can arise nonperturbatively, which we illustrate by considering effective string models. 
The nonperturbative contribution to the zero point energy vanishes in the Nambu model, but is 
nonzero when terms for extrinsic curvature are included. At one loop order, the nonperturbative 
contribution to the zero point energy is negative, regardless of the sign of the extrinsic curvature 
term. 

PACS numbers: 



I. INTRODUCTION 



The Wilson loop has a privileged status in a gauge theory. Although a nonlocal and 
composite operator, with dimensional regularization smooth loops in 3 + 1 dimensions are 
rendered finite by the usual renormalization constants for the gluon [1]. 

Lattice regularization, however, introduces an additional divergence. Considering the 
Wilson loop as the propagator for an infinitely heavy test quark, the new divergence, Eq, 
represents mass renormalization for the test quark. Eq has dimensions of mass, so if "a" is 
the lattice spacing, Eq ~ l/a, and multiplies the length of the loop. Extracting the quark 
antiquark potential from the Wilson loop in the usual manner, Eq corresponds to the zero 
point energy of the potential at asymptotically large distances. Although the (bare) quark 
antiquark potential has been measured with precision [2, 3, 4, 5], Eq is usually ignored. 

It is thus of interest to know how to renormalize Eq on the lattice. This is especially 
important for thermal Wilson loops, where in a pure gauge theory the Polyakov loop is the 
order parameter for deconfinement. While the bare Polyakov loop vanishes in the continuum 
limit, the renormalized loop does not [6, 7, 8, 9, 10, 11, 12, 13, 14]. However, as Eq varies, 
the the renormalized Polyakov loop I changes as I — ► e" 

-E /T t 

see, e.g., appendix C of 

ref. [10]. This nonperturbative ambiguity cannot be fixed by appealing to the perturbative 
regime at high temperature. 

In the first section of the paper, we follow previous analysis [6, 7, 9, 10] and suggest that 
there is a natural way to renormalize the Wilson loop such that the perturbative contribution 
to the zero point energy vanishes identically, EQ Srt = 0. 

This does not exclude nonperturbative contributions to the zero point energy. We illus- 
trate this by considering Eq 011 in one possible model of confinement, that of string models. 
If cr is the string tension, then just on dimensional grounds, Eq ou ~ y/a is possible. Even 
so, in the Nambu model [15, 16], or variants thereof [4, 17, 18], Eq ou = 0. 

A nonzero value of Eq 011 is generated only by models with massive modes on the world 
sheet. This arises by adding terms for the extrinsic curvature of the world sheet. A simple 
computation shows that at one loop order, Eq 011 is nonzero and negative, whether the cou- 
pling for the extrinsic curvature term is positive [19, 20, 21, 22] or negative [23, 24, 25, 26]. 

On the lattice, numerical simulations find that the quark antiquark potential appears to 
agree well with the simplest Nambu model down to rather short distances [3] . This suggests 
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that on the world sheet, any massive modes are heavy. In this case, the effective string 
theory is strongly coupled, so that the results of a one loop computation are only suggestive. 

II. RENORMALIZED ZERO POINT ENERGY 

Consider a rectangular Wilson loop of length t tot and width R. When t tot 3> R, the 
vacuum expectation value of the bare Wilson loop can be used to define the quark antiquark 
potential, 

(W) = (exp (ig j> A» dx^j J = exp (-V(R) ttot) • (1) 

At large distances, 

rv 

V(R^oo)~aR + E -- + ... (2) 

R 

Here a is the string tension, E is the zero point energy, and "a" is a constant. 

Assume that the theory has no dynamical quarks and confines, with the loop in the 
fundamental representation. Then o ^ 0, and a is universal, = 7r/12 in four spacetime 
dimensions [2]. In quantum mechanics, the value of the zero point energy simply produces 
a phase which multiplies the wave function, and is of no physical consequence. 

The potential above is a bare quantity, but only the zero point energy is ultraviolet 
divergent. The Wilson loop is related to the propagator for a massive test quark, with the 
zero point energy the additive shift in the mass. As the mass of the test field goes to infinity, 
the integral for this mass divergence is over three, instead of four, dimensions: 

The constants C\, c%, etc. depend upon the representation of the test particle, the details of 
lattice discretization, etc. 

The possibility of a zero point energy was recognized when the renormalization of the 
Wilson loop was first considered [1]. In perturbative computations it is of little consequence, 
because then dimensional regularization is natural, and the integral in Eq. (3) automatically 
vanishes, as a purely powerlike divergence in an odd number of dimensions. This remains 
true with other gauge invariant regulators, such as higher derivatives or Pauli-Villars, since 
they also eliminate such power law divergences. 
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To define a renormalized loop on the lattice, we introduce the renormalization constant 
Zti, which is a function of the (bare) coupling, g 2 : 

W£ arc = Z n {g 2 ) L / a W T Ti n . (4) 

With L is the length of the loop, L/a is the number of links for the path. 

One way of extracting the renormalization constant Z-ji(g 2 ) is to compute for different 
values of the lattice spacing, a, holding the physical length, L, fixed [7, 10]. This can be 
done directly from the numerical simulations, without using the perturbative expansion of 
Eq. (3). 

Another way of computing Z-n(g 2 ) is to compare the potential at short distances, V(R) 
as R — > 0, to the result in perturbation theory [6, 9, 10]. In a pure 577(3) gauge theory, the 
two methods agree with one another to the numerical accuracy tested [10]. 

Under this renormalization, it is possible to redefine 

Zn{9 2 ) - M9 2 ) L/a e-^ L /Zo , (5) 

so that the renormalized Wilson loop becomes 

WT -> e +S o eHL Z W£ n . (6) 

In these expressions Z Q is a pure number, which just shifts the overall normalization of 
the loop. We can always choose this normalization by considering very small loops; for 
Polyakov loops, this corresponds to the limit of very high temperature. For small loops, 
perturbative corrections are computable as a power series in ~ g 2 {L) ~ l/log(L). We can 
thus eliminate Zq by requiring that the small loops, suitably normalized, approach unity as 
L -> 0. 

Less trivial is the change due to Eq G , which shifts the renormalized loop by exp(i?Q ert L). 
This cannot be eliminated considering small loops, L — * 0, as represents a correction 
in a power of L, which in an asymptotically free theory is nonperturbative. 

We suggest that under perturbative renormalization, the only consistent choice is to take 
i?Q Crt = 0. We can renormalize the loop at zero temperature, where gluons are massless. For 
massless fields, the linearly divergent integral is uniformly proportional to the ultraviolet 
cutoff, which is ~ l/a. The basic point is that there are no finite terms ~ a . 
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This depends crucially upon the fact that at zero temperature, gluons are massless order 
by order in perturbation theory. Consider how the integral in Eq. (3) changes if the gluons 
did have a mass, m: 

ESrt™)--? j ^]^2=-9 2 (c 1 1 -- e 2 m + e 3 m 2 a + ...) , (7) 

for some constants ei, &3, etc. The ultraviolet divergent term, ~ 1/a, is independent of the 
mass m, and so the coefficient c% is the same as in Eq. (3). When m^O, though, there is 
a term ~ e2m which contributes a finite amount to the zero point energy. There are also 
terms at higher order which vanish as a — > 0, ~ e^m^a, etc. 

Our point is simply that when m — 0, then e 2 = 0, and there is no contribution to the 
zero point energy. This can be made more general. Consider first a non-Abelian gauge 
theory without dynamical quarks, and let the renormalization mass scale of the theory be 
A ren . To be invariant under the renormalization group, i?Q 6rt can only depend upon A rcn 
as EQ Brt ~ A rcn exp(— J dg/j3{g 2 )), where f3(g 2 ) is the /3-function for the theory. Doing so, 
however, means that E^ er is a dynamically generated mass scale. Such a nonperturbative 
mass scale is inconsistent with perturbative renormalization. 

Strictly speaking, this argument fails in the presence of dynamical quarks, where one could 
have Eq SV ~ m qua rk, with m quar k a current quark mass. We suggest that i?Q Crt vanishes even 
in the presence of dynamical quarks. However quarks modify the gluon propagator, the 
perturbative contribution to the zero point energy arises from a linearly divergent integral 
over a massless field, as it does in the pure glue theory. 

The bare zero point energy, Eq. (2), is then a sum of two contributions, 

E = El cn + £ non . (8) 

The perturbative contribution is ultraviolet divergent, with aEQ CTt ^ as a — > 0. This 
does not exclude nonperturbative contributions, Eq ou , for which aE™ 11 — > as a — > 0. We 
compute i?Q 0n in string models in the next section. Our point is simply that EQ ert can be 
uniquely determined: there is no freedom to vary i?Q ert by a finite amount by introducing 
£ pcrt . 

In particular, this implies that it is not possible to choose a E'q 61 * so that the renormalized 
potential, V(R) vanishes at a given distance [8]. Of course the renormalized potential will 
vanish at some distance Rq, but this cannot be fixed a priori. Instead, if follows from the 
renormalized quark antiquark potential. 
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There are several notable examples when the gluon develops a "mass" , and a finite zero 
point energy arises perturbatively. In all of these cases, however, renormalization proceeds 
as for m = 0; there is no ambiguity in how to compute such terms. 

The first is at nonzero temperature, where a Debye mass moebye ~ gT arises for Aq. 
Comparing to Eq. (7), this generates a zero point energy ~ gmoebye ~ g 3 T [11]. The 
coefficient €2 is positive, so at asymptotically high temperatures, the renormalized Polyakov 
loop approaches unity from above. A related phenomenon is familiar in Coulombic plasmas. 

Another case where a gauge field develops a mass is when it is coupled to a scalar field, 
<f), which then undergoes spontaneous symmetry breaking. Then the gauge fields acquire a 
Higgs mass, mHi ggs , and there is a zero point energy ~ g 2 mmggs- 

The last example where a perturbative zero point energy arises is for a gauge theory in 
three spacetime dimensions. Then the gauge coupling, g^ d , has dimensions of mass, and 
there can be a perturbative contribution to the zero point energy ~ g\ d . This cannot be 
distinguished from the nonperturbative contribution, since the square root of the string 
tension is also ~ g\ A . In this sense, it is much cleaner separating E^ OXi from E^ ert in four, 
instead of three, dimensions. 

In fact, in three dimensions one cannot avoid a zero point energy proportional to g 2 d 
[1, 7]. For a smooth Wilson loop in three dimensions, the only ultraviolet divergence is at 
one loop order; instead of Eq. (3), it is 

This then implies that bare loops vanish in the continuum limit as 

Zfigla) ~ exp (- Cue? g\ d L log (L/a)) ~ [jj . (10) 

As in four dimensions, in three dimensions bare loops vanish in the continuum limit, a — > 0. 
The suppression in three dimensions is only by a power of a, though, and is much weaker than 
the exponential suppression seen in four dimensions. This is seen in numerical simulations, 
where the bare Polyakov loop approaches unity at temperatures as low as several times the 
critical temperature [27]. 
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III. STRING MODELS 



To epitomize how a nonperturbative zero point energy can arise, in this section we con- 
sider a string model of the flux sheet. Of course there are many other models of confinement; 
presumably generic models also generate E^ on ^ 0. 

We consider a string action 

S = J d 2 z ^g ab (V a xJ (a + ^V 2 + ^V^ {V b x,) . (11) 

Here T> a are covariant derivatives with respect to the induced metric g a b = d a x^ on the 
surface x^(z); V 2 x' M = l/y/g d a (y/g g ab dbX^). 

In Eq. (11), u is the string tension, with dimensions of mass squared. The coupling for the 
extrinsic curvature term, k, is dimensionless. Lastly, the coupling A also has dimensions of 
mass squared. This action can be considered as the first terms in a power series of covariant 
derivatives on the world sheet. It will be clear later that our qualitative conclusions will not 
be altered by the presence of higher terms. 

A. Nambu model 

The simplest model is that for which k = A = oo, so the action only involves the area of 
the world sheet. In four dimensions, the exact solution for the potential is [16] 

V(R) = a Jr 2 -— « aR - — ^— + . . . . (12) 

v ; V 6(7 12R 288 a R 3 K J 

The term ~ 1/R is universal [2, 4, 20]. The potential is imaginary when R < ^n/(6a), 
which reflects the inconsistency of the pure Nambu model in four dimensions. 

For our purposes, all we need to recognize is that at large distances, the potential is aR 
times a power series in l/(aR 2 ). Thus there is no zero point energy in the Nambu model, 
E™ n = 0. 

We shall see in the next subsection that this arises because there are only massless modes 
in the Nambu model. Thus Eq ou remains zero for the models of Polchinski and Strassler 
[17] and of Liischer and Weisz [4]. While these models are rather different, in both higher 
derivative terms, like the extrinsic curvature term of Eq. (11), are treated as perturbations 
to the Nambu model. Thus there are only massless modes in such models, and Eq ou remains 
zero. 
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B. Rigid strings 



We next consider the model with positive sign for the coupling to the extrinsic curvature, 
k > 0, and neglect the coupling for a higher derivative term, A = oo. 

We compute for small fluctuations about a flat sheet, and follow Alvarez [15]. We let the 
number of spacetime dimensions, d, be arbitrary. At one loop order the d — 2 transverse 
fluctuations contribute to give an effective action, 

S eS = tr log (-<9 2 ) {-d 2 + m 2 ) ; -d 2 = -d 2 - d 2 , m 2 = k a . (13) 

The time is continuous, so the corresponding momenta uj is continuous. In the spatial 
direction, the transverse directions vanish at the end of the flux sheet, x^ r (0) = x^ T (R) = 0. 

Thus the associated momenta are discrete, p = nir / R, for n — 1, 2 

The effective Lagrangian involves a product of a massless mode, from the Nambu model, 
and a massive mode. The integral over the massless mode generates the term ~ 1/R. To 
do the integral over the massive mode, one can use analytic regularization [15], or consider 
the derivative with respect to m 2 : 

dS eS fd-2\ 1 fd-2\ A f duj 



tr — — = t, 



dm 2 \ 2 I —d 2 + m 2 \ 2 J ^ I 2n u 2 + (nir I R) 2 + m 2 

The uj integral can be done either directly, or by contour integration: 
dS eS fd-2\ A 1 



ft 



dm 2 ° V 4 J ^ (uj 2 + {mr/R) 2 + m 2 ) 1/2 ' 
Integrating with respect to m 2 , we obtain 

SMm > )=U J^)±((^)\,A'' 2 . (16) 

^ ' n=l ^ ' 

The sum over n is highly divergent, but can be done by using ^-function regularization. 

For m — 0, there is a single term, ~ 1/R. The sum over n gives £(— 1), and generates 
the Liischer term. 

When m ^ 0, we can expand the power series in powers of nit / R. The first term is 
proportional to m times ((0) = —1/2, so that the nonperturbative zero point energy is 

E r = - ( d ^ 2 ) ■ (17) 
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This result was obtained previously [21, 22], by computing in a large d expansion in 
weak coupling, k <C 1. Braaten, Pisarski, and Tse computed to leading order, and obtained 
Eq. (17), Eq. (16) of [21]. Braaten and Tse then computed at next to leading order in k, 
and so determined the corrections ~ Klog(m) to Eq. (17), Eq. (4.24) of [22]. 

Expanding S e g in powers of the spatial momenta, (nn/R) 2 , in principle we would expect 
corrections to the potential proportional to ~ 1/R 2 , ~ 1/R 4 , and so on. However, at 
one loop order a term ~ l/(R 2 ) e is proportional to £(— 2£); this vanishes when i is an 
intee "trivial" zero of the (^-function. Consequently, at one loop order there is only 

a contribution to the zero point energy, as all other corrections to the potential vanish. It 
is amusing to contrast this to the massless case, where simply on dimensional grounds the 
only contribution to Eq. (16) is ~ 1/R. 

That corrections ~ 1/(R 2 Y vanish is special to one loop order. This can be seen by the 
computations at next to leading order in k at large d, where from Eq. (4.23) of [22], at large 
distances the potential is 

ir .„ . d 7rd 1 7c 2 d 2 1 ii 2 d 2 1 

V(R -> oo) « o-R- -m + + ... (18) 

v ; 4 24 R 1152 o-R 3 384 maR* y 1 

The result for the Nambu model is obtained by taking d — > d — 2; then for d = 4, the terms 

~ 1/R and ~ 1/R 3 agree with Eq. (12). It is noteworthy that even at next to leading order 

in k, that there is no correction ~ 1/R 2 . We do not know if this is peculiar to next to 

leading order, or persists at higher order. 

There is a nonzero contribution ~ 1/R A , which can be understood as follows. At leading 

order, Eq. (16), there is a term ~ l/(m 3 i? 4 ) times £(—4), which vanishes. At next to leading 

order there is a contribution k times ~ l/(m 3 R A ), or ~ l/(maR 4 ), which is nonzero. 



C. Confining strings 

We turn next to the case of confining strings, [24, 25, 26]. We take this to mean a theory 
with negative coupling constant for the extrinsic curvature, k < 0, and A > 0. The effective 
Lagrangian is 

S cS = (^) tr log (-d 2 ) ((-d 2 ) 2 - 2M 2 (-d 2 ) + M 2 4 ) , (19) 

where 

2M 2 = - , M 2 4 = Aa . (20) 
9 



As usual, the massless mode generates the usual term ~ 1/R. 
The poles of the massive propagator are at 



p 2 = Ml ± iy/M* - Mf = Ml e ±2i * , tan(20) = J (j^ - 1 . (21) 

For the model to be physical, it is necessary that there are no poles on the real axis [25]. 
This gives the constraint, 

M%>Mf , * 2 >^. (22) 
As in the previous subsection we compute the derivative of the effective Lagrangian with 



ri 



respect to M 2 , 

dS cS fd-2\ 1 

1 tr -. — — T/ — ■ . . (23) 



<9M 2 4 \ 2 J (-d 2 ) 2 - 2MK-8 2 ) + M 2 4 
In taking this derivative we can assume that M\ is independent of M 2 . 

We first compute the zero point energy, neglecting the dependence upon the spatial 
momentum. This implies that the sum over the discrete spatial momenta generates ((0) = 
-1/2. 

It is then necessary to perform the integral over energies, u. For rigid strings, the propa- 
gator has two poles, at ±zm. In the present case there are four poles. Two are in the upper 
half plane, M 2 e 7,9 and M 2 e^ n ~ e \ and two in the lower half plane, M 2 e~ td and M 2 e~^ n+e \ In- 
cluding the contribution just of the two poles in the upper half plane, the pole at M 2 e l6 gives 
a residue ~ e~ j6, /(M| sin(26 1 )), while that at —M 2 e~ l6 gives a residue ~ e +ie /(M| sin(20)). 
Hence, 

dS eS . (d-2\ 1 



9M 2 4 tot V 16 J sin(^) M| ' 



Using sin(6 l ) = a/(1 — Mf /M 2 )/2, we find the nonperturbative contribution to the zero 
point energy to be 



= - (j^) \/M 2 - Ml . (25) 

This vanishes when M 2 = M%, but at this point the propagator has poles for real, Euclidean 
momenta, and the theory is not well defined [25]. 

Having obtained this result, one can also go through the algebra to determine corrections 
to the potential involving higher powers of 1/R 2 . When p = nir/R is included, the poles in 
the propagator are shifted. Even so, then only way that the spatial momentum p enters is as 
p 2 . Thus if one expands the result in powers of p, one will find that corrections are integral 



d 



10 



powers in p 2 . At one loop order, as for the rigid string this only involves the trivial zeroes 
of the ^-function, so that once again, the zero point energy is the only contribution to one 
loop order. As for the rigid string, we do not expect that this persists to higher loop order. 

Clearly one could consider higher derivative terms for an effective string model. By 
a similar analysis one expects that at one loop order, the only correction to the Nambu 
potential is E$ on < 0. 

IV. CONCLUSIONS 

We argued in Sec. II that a renormalized quark antiquark potential can be obtained from 
numerical simulations on the lattice, and that perturbative contributions to the associated 
zero point energy vanish. Nonperturbative contributions were computed in effective string 
models in Sec. III. They vanish in the Nambu model, but arise for either rigid or confining 
strings, with Eq ou < at one loop order. 

Numerical simulations appear to find that at large distances, corrections to the Nambu 
model are small [3]. This suggests that on the world sheet, any massive modes are heavy 
on the scale of yfo. This is a regime of strong coupling in effective string models, and so 
the one loop result for E™ a is not definitive. It will be interesting to see what numerical 
simulations find for both the sign and magnitude of Eq ou . 
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